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We estimate the possible variations of the gravitational constant G in the framework of a generalized (Bergmann- 
Wagoner-Nordtvedt) scalar-tensor theory of gravity on the basis of the field equations, without using their special 
solutions. Specific estimates are essentially related to the values of other cosmological parameters (the Hubble and 
acceleration parameters, the dark matter density etc.), but the values of G/G compatible with modern observations 
do not exceed 10~ 12 . 



1. Introduction 

Dirac's Large Numbers Hypothesis (LNH) is the ori- 
gin of many theoretical explorations of time- varying G . 
According to the LNH, the value of G/G should ap- 
proximately coincide with the Hubble rate. Although it 
has become clear in the recent decades that the Hubble 
rate is too high to be compatible with experiment, the 
enduring legacy of Dirac's bold stroke is the acceptance 
by modern theories of non-zero values of G/G as being 
potentially consistent with physical reality. 

There are three problems related to G , whose origin 
lies mainly in unified model predictions: 1) absolute G 
measurements, 2) possible time variations of G, 3) pos- 
sible range variations of G, i.e., non-Newtonian, or new 
interactions. For 1) and 3) see jl. 

After the original Dirac hypothesis some new con- 
cepts appeared and also some generalized theories of 
gravitation admitting variations of the effective gravita- 
tional coupling. We can single out three stages in the 
development of this field: 

1. Study of theories and hypotheses with variations 
of fundamental physical constants, their predic- 
tions and confrontation with experiments (1937- 
1977). 

2. Creation of theories admitting variations of an ef- 
fective gravitational constant in a particular sys- 
tem of units, analyses of experimental and ob- 
servational data within these theories jj] (1977- 
present). 
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3. Analyses of variations of fundamental physical 
constants within unified models Q (present). 

Different theoretical schemes lead to temporal vari- 
ations of the effective gravitational constant: 

1. Empirical models and theories of Dirac type, where 
G is simply replaced with G(t) . 

2. Numerous scalar-tensor theories (STT) of Jordan- 
Brans-Dicke type, with G depending on the scalar 
field <p(t) or a number of scalar fields. 

3. Gravitational theories with a nonminimally cou- 
pled (in particular, conformal) scalar field arising 
in different approaches Q (they can actually be 
treated as special cases of STT). 

4. Multidimensional unified theories in which there 
are dilatonic fields and effective scalar fields ap- 
pearing in our 4-dimensional spacetime from extra 
dimensions pj . They may also help one in solving 
the problem of a variable cosmological constant 
from Planckian to present values and the cosmic 
coincidence problem. 

A striking feature of the present status of theoretical 
physics is that there is no satisfactory theory unifying all 
four known interactions; most modern unification theo- 
ries do not admit unique and universal constant values 
of physical constants and of the Newtonian gravitational 
coupling constant G in particular. In this paper we dis- 
cuss the bounds that may be suggested by a general class 
of STT. One can mention that STT are among the vi- 
able alternatives to general relativity; on the one hand, 
they are widely used for comparison with observations 
and, on the other, their different versions emerge in the 
field limits of the candidate "theories of everything" . 
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Although the bounds on G and G(r) are in some 
classes of theories rather wide on purely theoretical 
grounds, since any theoretical model contains a number 
of adjustable parameters, we note that observational 
data concerning other phenomena, in particular, cos- 
mological data, place limits on the possible ranges of 
these adjustable parameters. 

Here we restrict ourselves to the problem of G (for 
G(r) see [1-4]). We show that various theories predict 
the value of G/G to be 10 _12 /yr or less. The signif- 
icance of this fact for experimental and observational 
determinations of the value of or upper bound on G 
is the following: any determination with error bounds 
significantly better than 10~ 12 /yr (combined with ex- 
perimental bounds on other parameters) will typically 
be compatible with only a small portion of existing the- 
oretical models and will therefore cast serious doubt on 
the viability of all other models. In short, a tight bound 
on G, in conjunction with other astrophysical observa- 
tions, will be a very effective "theory killer" and/or sig- 
nificantly reduce the class of viable theories. Any step 
forward in this direction will be of utmost significance. 

Some estimations of G had been done long ago in 
the framework of general scalar-tensor and multidimen- 
sional theories using the values of cosmological param- 
eters (Q, H, q etc) known at that time [Q, [|, ||, 9. 
It is easy to show that for modern values they predict 
G/G at the level of 10 _12 /yr and less (see also recent 
estimations of A. Miyazaki [Q, predicting time varia- 
tions of G at the level of 10~ 13 yr~ 1 for a Machian-type 
cosmological solution in the Brans-Dicke theory). 

The most reliable, by now, experimental bounds on 
G/G (spacecraft radar ranging ||) and lunar laser rang- 
ing ||) give a limit of 10~ 12 /yr), so any results at this 
level or better will be very important for solving the 
fundamental problem of variations of constants and for 
discriminating between viable unified theories. So, real- 
ization of such multipurpose new generation type space 
experiments like Satellite Energy Exchange (SEE) for 
measuring G and also absolute value of G and Yukawa 
type forces at the ranges of metres and the Earth radius 
PQ , become extremely topical. 

2. Scalar-tensor cosmology and 
variations of G 

We are going to estimate the order of magnitude of vari- 
ations of the gravitational constant G due to cosmologi- 
cal expansion in the framework of scalar-tensor theories 
(STT) of gravity, using modern data on the cosmological 
parameters. 

Consider the general (Bermann-Wagoner-Nordtvedt) 
class of STT where gravity is characterized by the metric 



g^ v and the scalar field cj> ; the action is 
S = J d 4 x^g{f(^n[g} 

+ h{<j ) )g^<j>, IJ ,ct>, v -2U{ct>) + L m }. (1) 

Here TZ[g] is the scalar curvature, g = \det(g^)\; /, h 
and U are certain functions of (p, varying from theory 
to theory, L rn is the matter Lagrangian. 

This formulation of the theory corresponds to the 
Jordan conformal frame, in which matter particles move 
along geodesies and hence the weak equivalence prin- 
ciple is valid, and non-gravitational fundamental con- 
stants do not change. In other words, this is the frame 
well describing the existing laboratory, geophysical and 
cosmological observations. 

Among the three functions of </> entering into ([!]) 
only two are independent since there is a freedom of 
transformations <p = (f>((j) new ) . We use this arbitrariness, 
choosing h(4>) = 1, as is done, e.g., in Ref. JTlJ . An- 
other standard parametrization is to put /(</>) = </> and 
h(<p) = Lo((f))/(j) (the Brans-Dicke parametrization of the 
general theory ([!])). In our parametrization h = 1 , the 
Brans-Dicke function lo(4>) is Ld(<f)) — f/f^; here and 
henceforth, the subscript 4> denotes a derivative with 
respect to <f>. The Brans-Dicke STT is the particular 
case u> — const , so that in (fil) 



/(</>) =</> 2 /(4w), h=l. 
The field equations that follow from ([l]) read 



(2) 
(3) 
(4) 



-<y^) + (v M v-^n)/- 

where □ is the D'Alembert operator, and the last term 
in (Q) is the energy-momentum tensor of matter. 

Consider now isotropic cosmological models with the 
standard FRW metric 

r 2 (d6 2 + sin 2 9 dep 2 )], 



ds 2 = dt 2 - a 2 {t) 



1 — kr 2 



(5) 

where a(t) is the scale factor of the Universe, and k — 
1, 0, —1 for closed, spatially flat and hyperbolic models, 
respectively. Accordingly, we assume 4> = <p(t) and the 
energy-momentum tensor of matter in the perfect fluid 
form T"^( m ) = diag(p, — p, — p, — p) (p is the density and 
p is the presume). 

The field equations in this case can be written as 
follows (the dot denotes d/dt: 



A 

a 



3 

-^(aa + a 2 + k) + Ua = 0, 



(6) 
(7) 



%(d 2 + k) = ±<p 2 + U-3*f + P , 
a z a 

4(2a5 + a 2 + k) = -U 2 +U - f -2-f -p. (8) 
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To connect these equations with observations, let us 
fix the cosmic time t at the present epoch (i.e., consider 
the instantaneous values of all quantities) and introduce 
the standard observables: 

H = a I a (the Hubble parameter), 
q = —ad/ a 2 (the deceleration parameter), 
— Pi Pa (the matter density parameter), 
where p cr is the critical density, or, in our model, the 
r.h.s. of Eq. @j in case k = : p CI = 3fH 2 . This 
is slightly different from the usual definition p cr = 
3H 2 /8irG where G is the Newtonian gravitational con- 
stant. The point is that the locally measured Newto- 
nian constant in STT differs from l/(87r/); provided 
the derivatives and f<p<p are sufficiently small, one 
has Ol 



87rG ff = 



1 2to 



(9) 



/2w + 3 

(more details can be found in Refs. |l4|, [lB) where the 
connection between G e s and uj was studied on the basis 
of cosmological solutions with local inhomogeneities and 
the equations of particle motion.) 

Since, according to the solar-system experiments, 
uj > 2500 , for our order-of-magnitude reasoning we can 
safely put 8irG = 1//, and, in particular, our definition 
of p CI now coincides with the standard one. 

The time variation of G, to a good approxiamtion, 

is 



G/G » -/// = gH, 



(10) 



where, for convenience, we have introduced the coeffi- 
cient g expressing G/G in terms of the Hubble param- 
eter H. 

Eqs. ([}])~(|^) contain too many arbitrary parameters 
for making a good estimate of g . Let us now intro- 
duce some restrictions according to the current state of 
observational cosmology: 

(i) k = (a spatially flat cosmological model, so that 
the total density of matter equals p C r); 

(ii) p — (the pressure of ordinary matter is negligible 
compared to the energy density); 

(iii) p = 0.3 p cr (the ordinary matter, including its dark 
component, contributes to only 0.3 of the critical den- 
sity; unusual matter, which is here represented by the 
scalar field, comprises the remaining 70 per cent). 

Then Eqs. (§) and (§) can be rewritten in the form 



U-3Hf = 2.1H 2 f, 



1 1? 



U-IH'f- f = (1 - 2q)H 2 f. 



(11) 
(12) 



Subtracting (|12|) from (|llj), we exclude the "cosmolog- 
ical constant" U, which can be quite large but whose 
precise value is hard to estimate. We obtain 



- Hf + / = (!.! + 2q)H 2 f. 



(13) 



The first term in Eq. ( |13[ ) can be represented in the 
form 

2 = f 2 (df/d^- 2 = /v/, 

and /// can be replaced with —gH . The term / can 
be neglected for our estimation purposes. To see this, let 
us use as an example the Brans-Dicke theory, in which 
/ = 4> 2 /{Auj) . We then have 

f + I <M\ 

here the first term is the same as the first term in 
Eq. (|l3|), times the small parameter 1/(2oj). Assuming 
that (jxf) is of the same order of magnitude as <f> 2 (or 
only slightly greater), we see that, generically, |/| <C <j> 2 ■ 
Note that our consideration is not restricted to the 
Brans-Dicke theory and concerns the model ([!]) with 
an arbitrary function /(</>) and an arbitrary potential 
U{<t>). 

Neglecting / , we see that ([l3]), divided by H 2 f , 
leads to a quadratic equation with respect to g : 



u>g 



9-q' = 0, 



(14) 



where q' = 1.1 + 2q. 

According to modern observations, the Universe is 
expanding with an acceleration, so that the parameter 
q is, roughly, —0.5 ± 0.2, hence we can take \q'\ < 0.4. 
(Note that this condition is only plausible rather than 
certain.) 

In case q' — we simply obtain g = —1/lj. Assum- 
ing 

H = h 100 ■ 100 km/(s.Mpc) w h wo ■ 10~ 10 yr" 1 

and uj > 2500, we come to the estimate 

\G/G\ <4-10- 1 % 00 yr-\ (15) 

where /iioo is, by modern views, close to 0.7. So (|l|) 
becomes 



\G/G\ <3-10^ 14 yr- 1 . 



(16) 



For nonzero values of q' , solving the quadratic equa- 
tion ([li]) and assuming q'u ^> 1 , we arrive at the esti- 
mate |<jr| ~ y/ q'/uj , so that, taking q' = 0.4 and again 
lj > 2500, we have instead of (fLa) 



\G/G\ < 1.3-10- 12 ^oo y^ 1 



0.9-10- 12 yr~\ 



(17) 



where we have again put hioo = 0.7. 

We conclude that, in the framework of the general 
STT, modern cosmological observations, taking into ac- 
count the solar-system data, restrict the possible vari- 
ation of G to values within 10~ 12 /yr. This estimate 
may be considerably tightened if the matter density pa- 
rameter Q m and the (negative) deceleration parameter 
q will be determined more precisely. 
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Our estimates are rather universal since they do not 
use special solutions to the field equations, but actually 
rest on the well-justified assumption that the expansion 
of the Universe occurs without abrupt changes in its 
parameters during a fairly long period before now. 

3. Discussion 

Summarizing the above considerations, we can conclude 
that restrictions on possible nonzero values of G give no 
bound on the possible class of generalized gravitation 
theories, but in the framework of some fixed theory any 
restriction on G restricts the possible class of models. 

We note that similar estimations of G can be made 
for different multicomponent multidimensional models 
giving a result on the level of 10~ 12 /yr and less 
for, e.g., dust and p-brane matter sources. 

We can also mention that the behaviour of the grav- 
itational constant can actually be much more complex 
and intriguing than a simple time (or even range) depen- 
dence: very recently, there appeared two papers, which 
may open a new series of theoretical and experimen- 
tal studies related to possible anisotropy in the absolute 
value of G (l(| and/or its possible dependence on the 
latitude and longitude of the laboratory where G was 
measured 17 . 



References 

[1] K.P. Staniukovich and V.N. Melnikov, "Hydrodynam- 
ics, Fields and Constants in the Theory of Gravity", 
Energoatomizdat, Moscow, 1983 (in Russian). 

[2] V.N. Melnikov, Int. J. Theor. Phys. 33, 1569 (1994). 

[3] V. de Sabbata, V.N. Melnikov and P.I. Pronin, Prog. 
Theor. Phys. 88, 623 (1992). 

[4] V.N. Melnikov, Gravity as a Key Problem of the Millen- 
nium. In: Proc. 2000 NASA/JPL Conference on Fun- 
damental Physics in Microgravity, CD-version, NASA 
Docu ment D-21522, 2001, p. 4.1-4.17, (Solvang, CA, 



USA); gr-qc/0007067 



V.N. Melnikov, Grav. & Cosmol. 6, 81 (2000). 

[5] K.A. Bronnikov, V.D. Ivashchuk and V.N. Melnikov, 
Nuovo Cim. B 102, 131 (1988). 

[6] N.A. Zaitsev and V.N. Melnikov, in: Problems of Gravi- 
tation and Elementary Particle Theory (PGEPT), 10th 
issue, Moscow, Atomizdat, 1979, p. 131; GR-9 Ab- 
stracts, Jena, 1980, p. 499. 

[7] A. Miyazaki, "Time Variation of the Gravitational Con- 
stant and the Machia n Solution in the Brans-Dicke The- 
ory", |gr-qc/0102003. 



[8] R. Hellings, Phys. Rev. Lett. 51„ 1609 (1983)] 

E.V. Pitjeva, in: "Dynamics and Astrometry of Natu- 
ral and Artificial Celestial Bodies", Kluwer Acad. Publ. 
Nethelands, 1997, p. 251. 

[9] J.O. Dickey et al., Science 265, 482 (1994). 



[10] 

[11] 

[12] 
[13] 

[14] 

[15] 

[16] 

[17] 



A. Sanders and W. Deeds, Phys. Rev. D 46, 480 (1992); 
V.N. Melnikov and A. Sanders, Ciencia Ergo Sum 8, 
357 (2001). 

B. Boisseau, G. Esposito-Farese, D. Polarski and A. A. 



Starobinsky, gr-qc/0001066 



A.J. Sanders et al., Clas. Qu. Grav. 17, 2331 (2000). 

V.N. Melnikov and V.D. Ivashchuk, in: Proc. JGR11, 
Tokyo, January 2002, Japan. 

V.N. Melnikov and A.G. Radynov, in: PGEPT, 14th 
issue, Moscow, Energoatomizdat, 1984, p. 73. 

V.N. Melnikov and A.G. Radynov, in: "On Relativity 
Theory", World Scientific, Singapore, 1985, p. 196. 

O.B. Karagioz et al., "Experimental Evidence that the 
Gravitational Constant Varies with Orientation" , Grav. 
& Cosmol. 8, No. 3 (31), 2002, to appear. 

J. P. Mbelek and M. Lachieze-Ray, "Possible Evidence 
from Laboratory Measure ments for Lati tude and Lon- 
gitude Dependence of G , gr-qc/0204064. 



